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Abstract 

In this paper, a new formulation for the three dimensional Euler equations is derived. 
Since the Euler system is hyperbolic-elliptic coupled in a subsonic region, so an effective 
decoupling of the hyperbolic and elliptic modes is essential for any development of the 
theory. The key idea in our formulation is to use the Bernoulli's law to reduce the dimen- 
sion of the velocity field by defining new variables (1, /32 — /^s — ^) and replacing u\ 

by the Bernoulli's function B through u? = txs2X«^- We find a conserved quantity for 
flows with a constant Bernoulli's function, which behaves like the scaled vorticity in the 
2-D case. More surprisingly, a system of new conservation laws can be derived, which 
is new even in the two dimensional case. We use this new formulation to construct a 
smooth subsonic Euler flow in a rectangular cylinder, which is also required to be ad- 
jacent to some special subsonic states. The same idea can be applied to obtain similar 
information for the 3-D incompressible Euler equations, the self-similar Euler equations, 
the steady Euler equations with damping, the steady Euler-Poisson equations and the 
steady Euler-Maxwell equations. 
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1 Introduction 

The three dimensional steady isentropic compressible Euler Equations expressing the 
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conservations of mass and momentum, take the following form: 



{pUl)x^ + {PU2)X2 + {P'^Zlxz = 0, 



{pul)xi + {pUlU2)x2 + {pUlU:i)a:^ + p. 







(1.1) 



< 



{pUlU2)xi + {pul)x2 + ipU2U3)x3 + P: 







{pUlUs)^;^ + (pU2ti3)a-2 + (puDxa +P: 







where u = {ui,U2,us) is the velocity field, p is the pressure and p is the density. We consider 
the ideal polytropic gas, which means that p{p) = Ap'^ , where ^ is a positive constant 
depending on the specific entropy and 7 G (1,3) is the adiabatic exponent. For simplicity, 

we take A = —. It follows from the momentum equations that the following important 

7 

Bernoulli's law holds: 



The original motivation to start such a research is the following: we try to construct a 
smooth subsonic flow in a three dimensional finitely long nozzle by formulating suitable or 
physically acceptable boundary conditions on upstream and downstream. Since the Euler 
system is hyperbolic-elliptic coupled in a subsonic region, the main task is to understand 
the coupling properties of the hyperbolic and elliptic modes in subsonic flows by finding a 
decomposition of the Euler system to effectively decouple the different modes. We consider 
this as a first step to get a close look at the difficult 3-D transonic shock problems in a 
divergent nozzle and hope that this may shed light on that problem. For the study of 
transonic shock in an infinitely long nozzle, one may refer to [HIT]. For the study of transonic 
shock problem in a duct or a De Laval nozzle, one may refer to O [U [lOl [TTl [131 1221 [Ml 
and [Ml [371 IHHl [39] for more details. This is a natural continuation of previous results in the 
two dimensional case discussed in [161 IMl [31] . 

There have been many relevant studies on the 2-D steady compressible Euler equations, 
especially for subsonic flows around a given profile and subsonic flows in an infinitely long 
nozzle. Profound understanding has been achieved both physically and mathematically for 
subsonic flows around a given profile by Frankl and Keldysh [20], L. Bers [H[5], M. Shiffman 
|28j and many other authors |141 \TE[ \T9\ \26\ I27j . By a variational method, Shiffman |28j 
proved that there is an optimal range [0, Mqo) for Mach number Mq at infinity that ensures 
the existence and uniqueness of a subsonic potential flow around a given profile with finite 
energy. Bers [1] presented a new proof of the existence and improved the uniqueness of a 
steady two-dimensional subsonic potential flow of a perfect gas around a given profile by 
removing the restriction of finite energy. In the recent works by Xie and Xin |[33l for 
subsonic potential flows in an infinitely long nozzle, they showed that there exists a critical 
value such that a global uniform subsonic flow exists uniquely in a general nozzle as long as 
the incoming mass flux is less than the critical value. They also obtained a class of subsonic- 
sonic flows by investigating the properties of these uniform subsonic flows and employing 
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the compensated compactness method. Such a strategy can not be apphed to the case in 
a multi-dimensional infinitely long nozzle (n > 3), but a different formulation involving the 
potential function has been employed by Du, Xin and Yan [T7] to obtain the existence and 
uniqueness results for subsonic potential flows in a multi-dimensional infinitely long nozzle 
(n > 2). The new ingredients of their analysis are methods of calculus of variations, the 
Moser iteration techniques for the potential equation and a blow-up argument for infinitely 
long nozzles. Concerning subsonic Euler flows in a 2-D infinitely long nozzle, Xie and Xin 
|35j explored the special structure of the 2-D Euler system and reduced the Euler system 
to a single second order quasi-linear elliptic equation, which yields the existence of subsonic 
flows when the variation of the Bernoulli's function is sufficiently small and the mass flux is 
in a suitable regime with an upper critical value. For subsonic Euler flow in an axisymmetric 
nozzle, one may refer to |15j . 

There are also some results doing this direction for both the 3-D compressible and incom- 
pressible Euler equations [TOl [HI EHl EHJ |39l HI [29] . In particular, Xin and Yin [39] developed 
a decomposition of the Euler equations for the uniqueness proof of the symmetric transonic 
shock solution and showed that the pressure satisfies a second order elliptic equation. They 
also studied the properties of the shock front. In [23j, they extended their techniques to 
obtain the uniqueness results in the case of variable end pressure. In [10], the author studied 
the transonic shock problems in 3-D nozzle with a general section, where he gave a decom- 
position of Euler equations in which the elliptic part and hyperbolic part is separated at the 
level of principle. To deal with the loss of derivative in the velocity field, the author first 
obtained the vorticity by solving transport equations, then the velocity field was obtained 
by solving an elliptic system with complex boundary conditions. However, it seems difficult 
to apply these techniques to construct a subsonic Euler flow in a 3-D bounded nozzle. In 
m , the author formulated an inflow-outflow problem for 3-D incompressible Euler system by 
prescribing the Bernoulli's function, the normal component of the vorticity field at the inlet 
of the flow region and the normal component of the velocity field on the whole boundary of 
the flow region. The velocity field is obtained by solving a first order elliptic system. Tang 
and Xin [29j have identified a class of additional boundary conditions for the vorticities and 
established the existence and uniqueness of solutions with non-vanishing vorticity for the 
three dimensional stationary incompressible Euler equations on simply connected bounded 
three dimensional domains with smooth boundary. 

In [16l [301 [31] ) we have characterized a set of physically acceptable boundary conditions 
that ensure the existence and uniqueness of a subsonic irrotational flow in a finitely long 
flat nozzle by prescribing the incoming flow angle and the Bernoulli's function at the inlet 
and the end pressure at the exit. In [16], we show that if the incoming flow is horizontal at 
the inlet and an appropriate pressure is prescribed at the exit, then there exist two positive 
constants niQ and mi with mg < mi, such that a global smooth subsonic irrotational flow 
exists uniquely in the nozzle, provided that the incoming mass flux m G [mo, mi). The 
boundary conditions we have imposed in the 2-D case have a natural extension in the 3-D 
case. One may also prescribe the incoming flow angles (— , — ) and the Bernoulli's function B 

Ul Ul 

at the inlet and the end pressure at the exit of the nozzle. We conjecture that such a problem 
should be well-posed for the 3-D Euler equations. These important clues help us a lot to find 
an effective decomposition for the 3-D Euler equations. It is natural to reformulate the Euler 
system in terms of new variables (s = lnp,j32 = —,(3-^ = —,B). Here one should note that 

Ul Ul 
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due to the Bernoulli's law ()1.2p . the Bernoulli's function B will possess the same regularity 
as the boundary data at the inlet, which is quite different from the velocity field. So the key 
idea in our new formulation in terms of (s, /32, /Ss, i?) is to use the Bernoulli's law to reduce 
the dimension of the velocity field by employing (l,/32 = — ,/?3 = — ) and replace ui by B 

Ul Ml 

through the simple algebraic formulae uf = — — ^2 ^^^ll ' ^^^^ '^^Y: we can explore the role 

of the Bernoulli's law in greater depth and hope that may simplify the Euler equations a little 
bit. In the case of no vacuum and stagnation points, the original Euler system is equivalent 
to (|2.2p . It is shown that s and 82/32 + d^P^ satisfy an elliptic system. At a first sight, one 
may regard (/32,/33) as hyperbolic modes since they satisfy transport equations (I4.10p . which 
will induce a loss of one derivative estimates when integrating along the particle path. Our 
first observation is that W = 82/33 — 83/32 + /3381/32 — /3281/33 is conserved along each particle 
path, at least for the case of flows with a constant Bernoulli's function. Indeed, W satisfies 
the equation (13. Gh . which behaves like the scaled vorticity in the 2-D case. And a physical 
interpretation of W is also given, see section [3] for more details. This, together with the first 
three equations in (12. 2j) may help to recover the loss of one derivative in the velocity field, so 
that the velocity field will possess the same regularity as the pressure. 

Our second observation is a system of new conservation laws, which is somehow surprising. 
The original motivation comes from our study on the 2-D Euler system in a curved nozzle 
(See [30j), where we use the flow angle = arctan — instead of the angular velocity w = —, 

Ml Ml 

to get a clear insight of the role played by the curvatures of the nozzle walls. We try to find 
a proper substitute in the 3-D case and choose the "sine function" of the flow angles. Finally 
a system of new conservation laws arises, whose physical interpretation is still unclear sofar. 
Even in the 2-D case, these conservation laws are new as far as we know. Yet, it seems that 
these information are still not enough to deal with the subsonic-sonic limit for a sequence 
of uniformly bounded 2-D Euler flows. How to use these new information to obtain some 
deeper understanding on the compressible Euler system is still under consideration. 

Back to our original problem, we can construct a smooth subsonic flow in a rectangular 
cylinder satisfying the prescribed incoming flow angle and the Bernoulli's function at the 
inlet and the given end pressure at the exit, which is also required to be adjacent to some 
special subsonic states. We emphasize that the background can have large vorticity. We 
develop an iteration scheme to obtain the solution as a fixed point of a contractive operator. 
The velocity field ^2 and /Sa will be solved by integrating along the particle path. The 
corner singularities are avoided by the reflection technique. The delicate nonlinear coupling 
between the hyperbolic modes and elliptic modes {/32,/33) makes it extremely difficult to 
have an effective iteration scheme in a 3-D nozzle with general section. How to develop an 
effective iteration scheme to obtain a smooth subsonic flow in a general 3-D nozzle is still 
under investigation. 

It is worthy noting that our idea can be applied to the 3-D steady incompressible Euler 
equations without any essential changes. That is to say, we can employ the Bernoulli's law to 
get a new decomposition for the 3-D incompressible Euler system. A new conserved quantity 
and a system of new conservation laws can be derived as for the compressible case. As we will 
see, the steady incompressible Euler system is always a hyperbolic-elliptic coupled system. 
We can construct a smooth incompressible Euler flow adjacent to some special states in a 
rectangular cylinder, which will satisfy the prescribed incoming flow angles, the Bernoulli's 



4 



function at the inlet and the end pressure at the exit. Our strategy can also be applied 
to the self-similar Euler equations, the steady Euler equations with damping, the steady 
Euler-Poisson equations and the steady Euler-Maxwell equations. 

The rest of the paper will be organized as follows. In the next section, we will present a 
new formulation for the compressible Euler equations and carry out a characteristic analysis 
to get an insight into the structure of the Euler system. The new conserved quantity and a 
system of new conservation laws will be derived in section [3l A smooth subsonic flow in a 
rectangular cylinder is constructed in section H] by employing the previous new formulation. 
In section [H we apply the same idea to give a new formulation of the incompressible Euler 
system. Some detailed calculations will be given in the appendix (section [6]) . 

2 A new formulation of the 3-D Isentropic Euler equations 



In this section, we will give a new formulation of the Euler equations and carry out a 
characteristic analysis for the new system to get some insights on the coupling structure of 
hyperbolic and elliptic modes in the compressible Euler equations. 

In the whole paper, it is assumed that the flow does not contain vacuum and any stagna- 
tion points, which means that p> and (^1,^2,^3) / (0,0,0) respectively. Without loss of 
generality, it is assumed that ui > 0. 



Define the following three new variables: /32 



U2 
Ul 



,/33 



U3 

Ul 



,s = In p. Multiplying the 



first equation in (11.11) by , dividing the z-th equation in (11. ID by puf and adding them 

puf 

together for i = 2, 3, 4, we obtain the following new system: 



dih + ^3^2 - hdzPz - hdis - (3ld2S + -^d2S - (32M3S = 0, 

uf 

dl^3 + (^282^3 - M2(32 - MlS - ^2^28 + -^dsS - ^lO^S = 0, 



-92/32 - 53/33 + 



c\p) 



dis - dis - 132823 - 133838 = 0, 



81B + hd2B + P383B = 
Using the third equation, we rewrite the above system as follows: 



(2.1) 



8is + I3282S + (33838 



M7 



818 + 82(32 + 83(33 = 0, 



8ih + 13282/32 + 13383/32 
8ih + (3282/33 + M3/33 



^/32a,. + ^92. = 0, 



—83S = 0, 



81B + (32d2B + P383B = 0. 



(2.2) 
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It is easy to verify that the above system ()2.2p is equivalent to (ll.ip . provided that the 
flow is smooth and does not contain vacuum and stagnation point. 

We focus mainly on subsonic flows in which the Euler equations become hyperbolic- 
elliptic coupled equations. Now we carry out a characteristic analysis to get an insight of the 
structure of the Euler equations. 

We rewrite (12.21) in a matrix form as 



Mi9iU + M2(92U + Mg^aU = 0. 



(2.3) 



Here U = [3,^2, 1^3, By and 

„2, 



Ml 







u 



,M2 



/32 1 



/33 1 



1 



Denote by A the root of det(AMi — ^2^2 — ^3-^3)- A simple calculation shows that 
det(AMi-6M2-e3M3) = (A-/3-e)'[(l— ^)A2-2A(/3-e)+(/3-e)' audit has 



/52 


1 










/33 





1 





c\p) 


h 













/33 








ul 








,M3 = 























h 











h 







ul 











/32 













P3 



one real root Xr = P ■ ^ with multiplicity 2 and two conjugate complex roots A^ = A/j it iXj, 



where Xr = — , A/ 

ai 



and ai = 1 ^"''^s = P • 4 = P26 + P343,a3 



ai 



Ul 



(/3 • if - ^fe^ + il). And aia3 - ai = ^[(^ - ml + Cl) " (/5 • 0^] > if the 

flow is subsonic, i.e. (?{p) > ul + u\ + u\. The corresponding left eigenvectors to Xr are the 
following 

1' = (0, 6 + /33 (/3 • , -6 - /32 (/? ■ , 0) , 1? = (0, 0, 0, 1) . 
The corresponding left eigenvectors to A^ are 



l^=lR± ^h = - a2, 6, 6, 0) ± ^{ " °^ 0, 0, 0)- 

ai 



The differential operators corresponding to Ir and 1/ are Di = (- 



-il32d2+f33d3), 82,83,0) 



'ul-c^{p) 

and D2 = (1, 0, 0, 0). The action of these two differential operations on ()2.2p yields an elliptic 
system for s and vj = 82(32 + 83(33: 
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dis + I32d2s + - -^dis + = 0, 



u 



1 ^ 

Ul — c [p) — c [p) 

3 3 (2.4) 

+ ^^-^TrT 5^(/3252 + /3353)A5iS + (52/32)' + 292/33^3/32 + (^3/33)' = 

Remark 2.1. T/ie differential operator corresponding to 1^ is D4 = (0,0,0,1). The action of 
D4 on ^2. S\) just yields the fourth equation in h2.S\) . indicating that the Bernoulli's function is 
conserved along the particle path. One may expect that the action of the differential operator 
corresponding to \\. would show us another conserved quantity for hyperbolic modes. Indeed, 
the differential operator corresponding to \\ is D3 = (0, 53 +/33(/3292 +/3353), — 92 — /32(/32C?2 + 
/3393),0) = (0,93-/339i+/33D,-92+/329i-/32D,0), where T> = 81+^282+^:^83. ApplymgTi^ 
to Ii2.2\) and tracing the principal term, one can find that Y)W is the principal term involved, 
where W = {P^di — 83)132 — {hdi — 82) (^3. However, we do not perform such a complex 
calculation here. It turns out that one just need to take the action o/(0, 83-/3381, — 92+/32<9i, 0) 
on \2.S^) . See the next section for more details. 

Remark 2.2. Here s and w = 82(32 + 83(33 satisfy an elliptic system, which is similar to the 
observation made in fW\. UT]/ . Indeed, the authors in fW\. 177] / carried out a characteristic 
analysis directly to the non-isentropic Euler equations and obtained a decomposition for the 
Euler equations. They defined a new variable v = 82U2 + 83U3 -82U1 ^i93Mi and showed 

Ul Ul 

that V and the pressure p satisfy the following elliptic system: 

81V + X1 {dl +8l)p + fi {Ui ,U2,U3,p,S)=Q, 

(2.5) 

V - X281P = f2iui,U2,U3,p,S). 

Here fi{ui,U2,U3,p, S) fori = 1,2 are functions of {ui,U2,U3,p, S) and the first and second 

V 

order derivatives of (ui,U2,U3,p, S). Simple calculations show that w = — . Our formulation 

Ul 

is different from 110[ liij / in the sense that we use the Bernoulli 's law to reduce the dimension 
of the velocity field by employing (l,/32 = — ,/33 = — ). Due to the Bernoulli's law 

Ul Ul 

the Bernoulli's function B will not lose regularity when integrating along the particle path, 
which is different from ui. We replace ui by B through the simple algebraic formulae uf = 



- h{p)) 
1 + /32 + /3I 



The 3-D Euler equations with a constant BernouUi's func- 
tion 
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3.1 A conserved quantity 



Before going into the details, we state the following simple but important fact. It is 
known that a two dimensional Euler flow is irrotational is equivalent to that the Bernoulli's 
function is a constant. However, this fact is not true in the three dimensional case. Indeed, 
by employing the following identity in vector calculus 

u • Vu = V(^|up) - u X cuHu, (3.1) 

and the conservation of momentum, we have 

V{^\u\^ + h{p)) -ux curlu = 0. (3.2) 

So if w = curlu = {wi,W2,W3) = 0, then V-B = and S is a constant. However, the 
constant Bernoulli's function only implies the vorticity field is parallel to the velocity field. 
One needs an additional condition such as d2U3 — d-iU2 = to guarantee that w = 0. In a 
word, the constant Bernoulli's function and ^2^3 — d-iU2 = will imply that the flow must be 
irrotational. However, d2U^ — d^U2 is not invariant under Galilean transformation. Also the 
three component of curlu are of the same status, there is no reasons to choose one of them 
specially. The following interesting calculation is due to the motivation that we try to find 
a better substitute of d2U^ — dsU2 in our new formulation. In this section, we concentrate 
on the flow with a constant Bernoulli's function, in which the vorticity field is parallel to the 
velocity field. 

A new observation is that the quantity W = 82/33 — 93/32 + /339i/32 — hdifi^ = (/Ss^i — 
93)132 — iP2di — 82)133 may serve as a good candidate. That is to say, = and the constant 
Bernoulli's function will imply that the flow must be irrotational. We first derive the equation 
satisfied by W and then give an interpretation of the physical meaning of W . 

Since the Bernoulli's function is a constant, thus = — ^ — To simplify the 

notation, we set G = 1 + /3| + /3| and define a new function K{s) of s satisfying K'{s) = 

— ; — — , then we rewrite the equations satisfied by Bo and /^^ as follows: 

liB-hie')) y HZ H6 

G-H8if32 + ^282^ + M^ih) - MiKis) + 82K{s) = 0, 



(3.3) 



G-'{di(33 + (3282^3 + Ms^s) - MiK{s) + 83K{s) = 0, 



Apply 13381 — 83 and — /32i9i + 82 to the above two equations respectively and add them 
together, one can show that W satisfies the following equation. 

BW + {82P2 + 83^3)W - G8iK{s)W - G'^WDG = 0. (3.4) 

It follows from the first equation in (j2.2p that W satisfies the following equation: 

nW - WUs - G-^WDG = 0. (3.5) 
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That is 

One may refer to the appendix for the detailed computations. 

Simple calculations show that = /3 • (V x /3) = — — - -\ — -— -| — -— = (u • w), 

Ul Ui Ul Ui Ml Ml uf 

here f5 = (l,/32,/33), which is the component of "scaled vorticity" along the "scaled velocity". 
Then ()3.6p becomes 

w 

The three components of curlu are functions dependent, so the "scaled helicity" — as a 

pG 

linear combination of wi,W2,W3, which is also invariant under Galilean transformation, seems 
to be a better candidate for the hyperbolic mode. 

Since the vorticity field is parallel to the velocity field, that is w = curlu = fj,(x)u, then 
we have 

uiwi U2W2 u^ws 2 , / xr- 

W = + + = p[x){l + + /^g) = ii{x)G. 

Ul Ul Ul Ul Ul Ul 

Then (|3.6p becomes 

d(^)=0. (3.8) 

Hence is conserved along the stream line, and plays a similar role as the "scaled vor- 

P 

(d u d Ul 

p 
0. 

Historically, the stationary solution of Euler equations with the vorticity field paralleling 
to the velocity field was called Beltrami flow and have been investigated for over a century. 
Arnol'd [2] had identified a class of flows with presumably chaotic streamline, that is the 
Arnol'd-Beltrami-Childress (ABC) flows with u = (j4 sin z+C cos y, B sm.x+A cos 2;, C sin y + 
Bcosx). Turbulent flows can be understood as a superposition of Beltrami flows has been 
proposed in the work of Constantin and Majda [12]. There are also some literatures on 
compressible Beltrami flows, for example [25]. Indeed, (13. 8p can be derived in the following 
simple way: using the divergent free condition for the vorticity field, curlu = fj,{x)u will 
imply that div{fi{x)u) = 0, then by employing the density equation, we directly obtain (13. Sp . 

Now we explain the reason why we spend so much time and energy on the above com- 
plicated calculations for (13. Sp . The point is our calculations also work for the isentropic 
Euler equations. For the isentropic Euler system, similar calculations as above show that W 
satisfies the following equation: 



BW -WDs-G-^WDG ^^^^^^ 



2{B - h{p)f 
+ hidzBdis - diBdss) + P^{diBd2S - d2Bdis) 



{d2Bdss - d^Bd2s) 
0. 



(3.9) 



That is 
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Remark 3.1. Our original goal is to construct a smooth subsonic flow satisfying the prescribed 
incoming flow angles and the Bernoulli 's function at the inlet and the given end pressure at the 
exit. Assume that the boundary data are in C^'". Since the Bernoulli's function is conserved 
and will not lose regularity when integrating along the particle path, then B G C^'". It follows 
from Ii2. 2\) that s satisfies a second order elliptic equation 

di - l)5is - f32d2S - M^s^ + 92 ( - P2dis + - f3l)d2S - Msd^s 

+ 53 - Mis - hP^d2S + - Pl)d^s^ - (^(^ - l)dis - hd2S - M^s^ 

{82(^2? + {ds^sf + 2^2/3353/35 

(3.11) 

This shows that s also possesses C^'" regularity in general. However, (^2 and (3^ will lose 
one order derivative in the process of integrating the transport equations. Roughly speaking, 
the quantity W may help to raise the regularity of (32 and 13^. Since W satisfies \3.1 W 
should have C^'" regularity in general. Combining this with the first three equations in li2.S\) . 
one can obtain a divergent-curl system for (32 and f3s, which may help to raise one more 
regularity for /32 and (3-^. So the velocity field would possess the same regularity as that of 
the pressure. The above informal argument indicates that whether the Bernoulli's function is 
constant or not is not so essential for enhancing the regularity of the velocity field. 

Remark 3.2. Indeed, one can also derive the equations \3. 6\) and ( [3. 7| ) directly from the 
original Euler equations It follows from the original Euler equations that the vorticity 

w satisfies the following equations: 

uidi'w + U2d2^ + U393W + divu • w — (w • V)u = 0. (3.12) 

Combine these with the density equation, we have 

{uidi + U2d2 + usds) - (^^) u = 0- (3-13) 

It follows from \3.1S\) that u • w satisfies 

(u • V)(u • w) + divvL • (u • w) - (w • V)||uf = 0. (3.14) 
Employing the Bernoulli's law and the density equation, one can derive the following equation: 

||u||2(u.V)(^)-^(u.V)(||uf) + ^^(||uf(w.V)p-(u.w)(u.V)p) =0. (3.15) 



That is 

^p||u||^y p^iju"" 



" • V) ( ^) + ( ||uf (w . V)p - (u • w)(u . V)p ) = 0. (3.16) 



IfVB = 0, that is u || w, then the last term drops out. We obtain (3.1) again. It should 
be emphasized that \3.15\) holds also for flows with stagnation points. Hence i/ u • w = at 
the inlet, it holds also in the whole nozzle. 
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Remark 3.3. We also remark here that 13. 2^ implies that the Bernoulli's function is conserved 
not only along the stream line, hut also the vortex line. The last conclusion seems not getting 
enough attention before. Indeed, since the vorticity field is divergent free, we obtain the 
following new conservation law: 

div{Bcurlu) = 0. (3.17) 



3.2 A system of new conservation laws 



In the following, we assume that the Bernoulli's function is a constant and the speed |u| 
has a positive lower bound 5. It is easy to see that G satisfies the following Riccati-type 
equation: 

Define three new variables: Ki = G '2,Ki = — ^ = PiKi,i = 2,3. Here K2 and 

G2 

behave like the "sine function" of "the flow angles" . In the following, we derive the equations 
satisfied by K2 and K3. To simplify the notation, we define two new functions I{p) = 

gJo 2t(B-h(t)) Qj^d Q{p) = / —77^ dt. Since the speed |u| has a lower bound 6, 

Jo 2(i^ - h{t)) 

f c^it) 1 

/ — ; , , , . < -p^hip), so Up) and Qip) are well-defined and Qip) > for p > 0. One 

Jq 2t{B - h[t)) 

^, , . „ I'(p) c^ip) 
should note that I[p) satisfies — 



i{p) 2p{B-h{p)y 

It follows from (j3.18p that Ki satisfies: 

While the first equation in (12. 2p yields 

—^P - TTTT^^TT^^r'^i/' + ^i(^2/32 + ^3/33) = 0. (3.20) 
p 2p{B - h{p)) 

Combing these two equations, we obtain 

:^Dp + DKi - -^M-— KiDp + i^i (92/32 + 93/33) = 0. (3.21) 
p 2p{B - h{p)) 

That is 

-K • Vp + div¥. - K ■ Vp = 0. (3.22) 

p 2p{B - h{p)) 

Using the definition of /(p), we have the following new conservation law, which is similar to 
the conservation of mass: 

div{pA) = 0. (3.23) 

where A = -i-K = -l-(Ki, K2, K^). 

Hp) Hp) 
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Similarly, one can also derive the equations satisfied by i^j, i = 1, 2, 3, as 



Di^3 



cHp) 



2p{B - h{p)) 

cHp) 

2p{B - h{p)) 
cHp) 



Ki-Dp + 
K2BP + 



c\p) 



2p{B - h{p)) 

cHp) 
2p{B - h{p)) 

cHp) 



That is: 



2p{B-h{p)) ^ 2p{B-h{p)) 
pK-VI{p) -^'-^ 



pK 




-Ki 


< pK 




-K2 


p-K 


VK3 





Hp) 


+ 


2{B - h{p)) 


pK ■ VI{p) 




cHp) 


i{p) 


+ 


2{B-h{p)) 


p-K ■ Vlip) 


+ 


cHp) 



K^'dip = 0. 

K^^d2P = 0, 
K^^dsp = 0. 

dip = 0. 
d2P = 0, 
dsP = 0. 



(3.24) 



(3.25) 



I{p) 2{B - h{p)) 

With the help of (j3.23p . we obtain the following conservation laws, which are similar to 
the conservation of momentum: 



di{pAl) + ^2(^1^2) + d^ipAiA^) + diQ{p) = 0, 
di{pAiA2) + d2{pAl) + 93(^2^3) + d2Q{p) = 0, 
diipAiA^) + d2{pA2A^) + d^{pAl) + dsQip) = 0. 
For the isentropic Euler equations, since (5i + /32d2 + /S^d^jB = 0, we still have 



(3.26) 



• Vp + divK - . ^ K ■ VI{p, B) = 0, 
P I{P,B) 



(3.27) 



rf . 

where I{p, B) = e''° MB-Kt)) Hence one can still obtain the following conservation law: 



rdt 



di P 



I{p,B] 



+ d2{p 



K2 



+ dAp 



I{P,B) 



0. 



(3.28) 



IiP,B), 

However, similar conservation laws as ()3.26p are not true in general. 

For the non-isentropic Euler equations, one can perform the same operations to get similar 
information. Actually, the non-isentropic Euler equations can be rewritten as 

-{dip + i32d2P + Mzp) - -\diP + d2h + 93/33 = 0, 

P puf 

dih + /3292/32 + - -\P2d1P + -\d2P = 0, 

pui puf 

dips + /3292/33 + M^iPs - -\MiP + -^dsP = 0, 

puf puf 

diB + /32d2B + p^d^B = 0, 
diS + f32d2S + /SsdsS = 0. 



(3.29) 
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1 1 
In this case, B = -{uf + + uj) + h{p, S) = + /3| + /3|) + S"), where h{p,S) 

dpPir, S) _ 2iB-hip,S)) 



r-E^dr. Hence.? 

Jo 'T 



10 I ^-t- P3 

Similar calculations as previous sections show that W = 82/33 — 83/32 + /339i/32 — /32'9i/33 
satisfies the following equation: 



D( ^)+ir-T77T7;V71^-T77TT^<iB-hiP,S))8spiP,S)P-iVPxVS) 



pGj 2pHP,S){B-h{P,S))\ 

+ p{P, S) [/3 • (VP X VB) - dsh{P, S)f3 ■ (VP x VS)] | = 0, 

where we have rewritten p and h as smooth functions of (P, S) . 
Also a similar conservation law holds: 

rp apP(t,s) . 
where I{p,B,S) = e-'o 2*(s-Mt.s))''*. 

Remark 3.4. Since G satisfies a Riccati-type equation icl.l8\) . one may expect some blow-up 
results. However, the ambiguous sign ofdis in the coefficient ofG^ in i3.18\) makes the whole 
argument nontrivial. One should note that G blows up means that the fluid will turn around 
in the flow region. 

Remark 3.5. The physical interpretations of i3. 23\) and 113. 26\) are not clear at this moment. 
How to use these new information to obtain some useful results is an interesting problem and 
will be investigated further. 

Remark 3.6. One may perform the same operation as we have done in Remark \3.3l to i3.23\) 
and \3.26\) to obtain a new conservation law involving curl A. However, one should note that 
the Bernoulli 's function for i3.26\) is just a constant, so no new information can be obtained 
in this case. 

Remark 3.7. One may perform similar operations as we have done for the original Euler 
equations, to the new conservation laws \3. 23\) and i3.26\) . it turns out that one will get 
f3.23\) and 1^3. 26\) again. That is, 13.23\) and /i3.26\) is invariant under our operations. 

Remark 3.8. These calculations certainly work for the two dimensional Euler system. Indeed, 
W will always be zero. We can obtain a system of new conservation laws as \3.23\} and f3.26\) 
for potential flows and a new conservation law as 113. 28\} or i3. 31\) for the 2-D Euler equations. 
Even in this case, all these information are completely new as far as we know. However, it 
seems that these new information are still not enough to deal with the subsonic-sonic limit 
for the 2-D Euler flow. 

Remark 3.9. We may apply the same idea to the self-similar compressible Euler equations. 
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which take the following form: 



(3.32) 



(pf^i)Ci + (Pf^2)6 + (pUs)^, +3/9 = 0, 

pUid^.Ui + pU2d^,Ui + pC/3%f/i + %p + pUi = 0, 

/)C/l%C/2 + />C/2%C/2 + + %p + pU2 = 0, 



i/eree = (^1,6,6) = (y'^.y) anrf (C/i, C/2, C/3)(e) = ^2, ^Z3)(e) - (6, 6, Cs)- T/ie 



mam difference in \3. 32\) is that the Bernoulli's law does not hold any more. One can still 
define 
for (EI 



define s = In p, (^2 = 77-5/33 = 7^, B = —{Uf + [/| + C/|) + h{p) and obtain a new formulation 
U\ U\ 2 



(%s + /32%s + /33ag3s) 



c\p) 



^2 -d^^s + 82^2 + 53/33 + ^ = 0> 



%/32 + /32%/32 + /33%/32 - ^/32%5 + = 0, 

%/33 + /32%/33 + /33%/33 - -T^m.s + = 0, 



(3.33) 



d^,B + hdi-^B + l33d(..^B + 



0. 



Although the Bernoulli's law does not hold in this case, the Bernoulli's function B still has 
the same regularity as that of the pressure, and W = 82/33 — 83/32 + /3381/32 — /32(?i/33 may also 
help to raise the regularity of /32 and /33. Hence the reformulation 113. 33\) may provide some 
new insights on the self-similar Euler equations. Same arguments work for the steady Euler 
equations with damping. 

Remark 3.10. In ]32^ . we will employ the same idea to investigate the structural stability 
of some steady subsonic solutions for 1-D Euler-Poisson system under multi- dimensional 
perturbations of suitable boundary conditions and background charge. 



4 Subsonic Euler flows in a rectangular cylinder 



In this section, we try to construct a subsonic Euler flow in a rectangular cylinder by 
imposing suitable boundary conditions at the inlet and exit. The cylinder will be = 
[0, 1] X [0, 1] X [0, 1]. Actually the flow we construct below will be adjacent to some special 
subsonic states {pq,uq{x2,X3),^,{)). Here po is a constant, uq{x2,X3) G C^'°([0, 1] x [0,1]) 

satisfying < ^0(2:2,^:3) < c{pq). We also set Bq{x2,X3) = -ul{x2,X3) + h{po),so = InpQ. 



14 



At the inlet of the nozzle xi = 0, we impose the flow angles and the Bernoulli's function: 

ft(0,X2,X3) = e/3r(x2,X3),i = 2,3, 

Here the compatibility conditions should be satisfied: 



(4.1) 



a^/3r(0,X3) = 9^/3^(1, 2^3) = 0, 
9^/3^(^2,0) = a|/3r(x2,l) = 0, j = 0,2, 

a2'i?*"(0,X3) = a|i?^"(l,X3) = afi?^"(x2,0) =ati?^"(x2,l) = 0, A: = 1,3. 
At the exit of the nozzle xi = 1, we prescribe the end pressure: 



(4.2) 



p(l,X2,X3) 



o7(so+ese(x2,a;3)) 
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(4.3) 



(4.4) 



(4.5) 



Here we also require that Se satisfies the following compatibility conditions: 

^Se(0,X3) = ^Se(l,X3) = 0, 

9^Se(x2,0) = 9^Se(a;2, 1) = 0, j = l,3. 
While on the nozzle walls, the usual slip boundary condition is imposed: 

li2(xi, 0,2:3) = ^2(2:1, 1,X3) = 0, 
li3(xi,X2,0) = U-i(xx,X2, 1) = 0. 



Mathematically, we are going to prove that (jl.ip with boundary conditions (|4.1|) . 
and (|4.5p satisfying compatibility conditions (|4.2p and (|4.4p . has a unique subsonic solution. 

4.1 Extension to the domain ^7e = [0,1] x 

Suppose that the flow {p,ui,U2,U3) € C'^'"(J^) x C^'"(d)^ to be constructed has the 
following properties: 



dl{p,ui){xi,0,X3) = dl{p,ui){xi, 1,X3) = 0, 
di{p,ui){xi,X2,0) = dl{p,ui){xi,X2, 1) =0, j = 1,3, 

(4.6) 

52M2(a:i,0,X3) = 92ti2(2;i,l,X3) = 0, 

dlu3{xi,X2,0) = dlu3{xuX2,l) =0, A: = 0,2. 
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Then we may extend (p, ui, n2, ua) in the fohowing way, to (/5, ui,U2, U3) € C'^'"([0, 1] x 
C2''^([0,l] X M2)3. 

For (xi,X2,X3) G [0,1] x [—1,1] x [—1,1], we define {p,ui,U2,U3) as fohows 



(p, Ui,n2,?i3)(x) 



(p,til,'U2,ti3)(xi,X2,X3), if (x2,a;3) G [0, 1] X [0, 1], 
{p,Ui, -U2,U3){xi, -X2,X3), if (X2,X3) G [-1,0] X [0, 1], 
{p,Ui,U2, -U3){xi,X2, -X3), if (X2,X3) G [0, 1] X [-1,0], 
{p,Ui,-U2,-U3){xi,-X2,-X3), if (X2,X3) G [-1,0] X [-1,0]. 



Then we extend {p,ui,U2,U3) periodicahy to [0,1] x with period 2. It is easy to verify 
that {p,ui,U2,U3) belong to C^'"([0, 1] x M^) x C^'"([0, 1] x M^)^. Moreover, {p,ui,U2,U3) 
satisfy the Compressible Euler equations. Due to these reasons, one may directly work on the 
domain [0, 1] x T^ (Here T^ is a 2-torus), so that the difficulty caused by corner singularities 
and slip boundary conditions can be avoided. We may extend (uq, Bq, pi^ , P"^^ , i?*", Se) to T^, 
which win stih be denoted by {uo,Bo, p^ , pf ^B'"" , Se). 

4.2 Main results 

The main result in this section is the following existence and uniqueness theorem. 

Theorem 4.1. Given 5*", Se) G C^'"(t2) satisfying the compatibility conditions 

(j4.2p and (j4.4p . there exists a positive small number eo, which depends on the background 
subsonic state (po, uo(x2, X3), 0, 0) and (/J™, Z?™, i?*", Se), such that if < e < eo, then there 
exists a unique smooth subsonic fiow (p, ui, U2, ?/3)g C^'"(r2e) x (C^'"(r2e))^ to (jl.ip satisfying 
boundary conditions ()4.ip . (j4.3p and (j4.5p . Moreover, the following estimate holds: 

\\{ui,U2,U3) - (uo(x2,X3),0,0)|[c2,«(Q^) + ||/0 - Po ||c3." (f^e) ^ C'^. (4.7) 

Here C is a constant depending on (po, no(a;2, xs), 0, 0) and {13^ , PI"' , B^"" , Se)- 

Since we are looking for a smooth subsonic Euler fiow in the whole nozzle, which is also 
close to some special subsonic states, no vacuum and stagnation points will appear in the 
fiow region, we may use the new formulation developed in the section [2] and section [3j Hence 
to prove the above theorem, it is equivalent to show that the following system has a unique 
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smooth subsonic solution: 



dis + P2d2S + ^3^35 - -^d^s + 82^2 + 03/33 = 0, 
diP2 + ^282^2 + PM2 - -^hdis + -^d2S = 0, 



dlk + ^2^2/33 + /3353/33 - ^^3^15 + ^^3^ = 0, 

ui ui 



cHp) 



(4.8) 



We can derive that s = s 
respectively: 



diB + ^32^25 + ^3^3^ = 0, 
(3i{0,X2,X3) = ePi"'{x2,X3), for i = 2,3, 
B{0,X2,X3) = Bq + eB'''^{x2,X3), and ni(0, a;2, X3) > 0, 

s(l, X2, X3) = So + ^Se{x2, X3). 

50? 02,^3 and B = B — Bo satisfy the following equations, 



di 



- l)diS - P2d2S - MSS^ + ^2 ( - p2dlS + (^ - Pi)d' 



+^3 ( - /335l^ - P2P3d2S + (^ - Pl)d3S 



U, 



c^(p) ' ^ 

(-^-l)aiS-/3292S--/3353^ 



+ (^2/32)' + {dsPsf + 252^3^3^2 = 0, 



(^ - 1)^15 - ^2^25 - P3d3S^ (0, X2, 2:3) = e(52^r + ^3/3^), 
s(l,a;2,X3) = €Se{x2,X3). 



(4.9) 



5l^2 + P2d2P2 + /9353/S2 - ^/325lS + ^d2S = 0, 



5l^3 + /3252/33 + P3d3l33 - ^PsdiS + ^838 = 0, 

uf uf 

/32(0,X2,X3) = e^^"(x2,X3), 
/33(0,X2,X3) = e^|"(x2,X3), 

V 

5i5 + ^32^25 + MsB = -l32d2Bo - MzBq, 

5(0,X2,X3) = e5^"(x2,X3). 



c\p) 



(4.10) 



(4.11) 
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4.3 Proof of Theorem 14.11 



The main idea is simple: we construct an operator A on a suitable space, which will be 
bounded in a high order norm and contraction in a low order norm. 
The solution class will be given by 

S = {{s,^2,^3,B) : ||s||c3,a(ne) < ll/32,/33,-B||c2,a(n^) < 6.} 

Here 6 will be determined later. For any given (s, (^2, (S^, B) G H, we define an operator 
A : {s, /32, (3^, B) i— > (s, (32, fH^, B) mapping from H to itself, through the following iteration 
scheme. 

Step 1. To obtain s by solving the following linearized elliptic system: 

^1 ( - l)dis - hd2S - hd^i] + ^2 f - hdis + - ~Pl)d2-s - hhd^s 



w 



+dA - hdis - hhd2S + { 



l)dis - hd2S - Mss (0, X2,xs) = e{d2PT + ^s/^D, 



Here Fi 



S(1,X2,X3) = eSe(x2,X3). 

c\p) 



(4.12) 



l)ail - - - (92/32)" + {ds/SsY + 202/33^3^2 , P = 6 



pS + So 



C^'"(r2) with the following estimate: 



1 + /3| + /3| ■ 

By the standard elliptic estimates, we can show that (|4.12|) has a unique solution s G 



< C2{S^ + e). 



(4.13) 



Here Ci depends only on the background solution and J^e, and C2 depends also on {/3^, /3i^, i?*", Se). 
Step 2. To obtain /32;/33 by solving the following hyperbolic equations: 

dlp2 + I32d2(32 + /3353/32 - -^Mis + -^d2S = 0, 



dik + /3252/33 + /33a3/33 - -^kdis + -^d3S = 0, 

/32iO,X2,X3) = e/3^"(x2,X3), 
/33(0,X2,X3) = e^f (X2,X3). 



(4.14) 
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The particle path (r, X2(7"; x), X3(t; x)) through (xi,X2,a;3), is defined by the following 
ordinary differential equations: 



(ix2(T;x) 
(ix3(r;x) 



^2(r, X2(r;x),X3(r;x)), 
P3{t, X2(r;x),X3(T;x)), 



(4.15) 



dr 

X2(xi;x) = X2, 

X2(xi;x) = X3. 

Since (/32,/33)GC^'"(^^e),it is easy to prove that (x2(0; x), X3(0; x)) belong to C^'"(Oe) with 
respect to x. Indeed, we have the following estimates: 

||(x2(0;x),X3(0;x))|b2..(f,^) <C73. (4.16) 

Then by the characteristic method, it holds that: 

/32(x) = e/3^"(x2(0; x), f 3(0; x)) + J^' (^^{^^dis - 823)^ (r, 52(t; x), X3(r; x))(ir, 

/33(x) = e/3r(x2(0; x), X3(0; x)) + (^^{^^dis - d^s)^ (r, X2(t; x), X3(r; x))(ir. 

(4.17) 

These enable one to obtain the following estimate: 

||(/32(x),/33(x))|b2,.(^^) < C74[e||(/3r,/3r)lb2,.(T2) + \\s\\c^,<^in)] < C,{6' + e). (4.18) 
Step 3. To resolve B. Solving 



diB + hd2B + MzB = -P2d2Bo - MzBo, 

5(0,X2,X3) =eB*"(x2,X3). 

As above, B has the following formulae 



(4.19) 



5 = eB*"(x2(0;x),X3(0;x)) + 



/3292-Bo - /SsdsBo (r, X2(t; x), X3(t; x))(ir. 



Furthermore, the following estimate holds: 

<e||i?*"||c2,<.(T2) + C5||(/32(x),/33(x))|lc2,.(n.) <C6(<5' + 6). (4.20) 
Collecting ah the estimates (jH3|), and ()i:20D gives 

I|s||c3,"(c,) <C7(52 + 6),||/32,/33,5|b2,.(f7.) <C7(<52 + e). (4.21) 

Here Cj = max{C2, C5, Cg}, which depends only on the background solution, fig and (/J^i /^s", -B* 
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Choose ei small enough, such that if e < ei, then C|e < ^. Set 6 = 20^6, then C^d < ^ 

and C'7{S^ + e) < Cye = 5. This implies that A maps S to itself. 

It remains to show that the mapping A : S ^ H is a contraction operator. Suppose 
A : ^ (s^/32^/3t,B^) for k = 1,2. Define the difference (Fi, Fa, ^3, n) = 

is' - - Pi Pi - Pi B' - Bl and (Yi, Y^, Y^) = (1^ - 1^, ~PI - Pi P^ - PI B' - B^). 

Step 1. Estimate of Yl. 

We rewrite the equation satisfied by s as follows: 
Y,di(aij{^,~P2,h.^)dA = (5Zai,(l,^2,/33,Ml) -({d2~P2?+{dzh?+'^d2hdz~pX 

i,j=l ^ ^ ^j=l ^ ^ ^ 

(4.22) 



Then Yi satisfies the following elliptic system: 



3 



i=l 



-Y,[aij{f,PlPl,B')djYi]{0,X2,xs) = Gi, 

yi(l,X2,X3) =0. 



Here 



Gi = J2[ Ki^'JlPhB') - aiji^'', PlPlB'')]djS^ , i = 1, 2, 3; 



1=1 



i=l 



H2 = {d2~Plf + {diPlf + 252^1^34' - Wl? + {d^~Plf + 2d2~Pld^~Pl 



Thus the following estimate holds: 

/ 3 

Il^l||c2.«(ae) <CA XI + ll(-^l'-^2)||c-(ne) + 1 1 1 1 CL" (T2) 

^ i=l 

< C95(^||yi||c2,.(o,) + IIY2, Y3,y4||ci-(ne) 
Step 2. Estimate of F2, Y3. 



(4.23) 



(4.24) 



(4.25) 
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It follows from (14.140 . I2 and satisfy the following system: 
51^3 + ^2d2Y3 + ^IdsY^i 

Y2{0,X2,Xs)=0, 
>3(0,X2,X3) =0. 



-^{Pld,-d2)Y,=K,, 
pP^{^ld,-d,)Y,=K2. 



(4.26) 



Here 



K2 



-(m/3| + y393/3|) + 



/3i 



02S , 



Then it follows directly that: 



\\{Y2,Y3)\\c^.^i^^^)<CiJ5\\{Yi,Y2,Y3,^^^^^ 



(He) 



< Cn,5 ||yi||c2..(ne) + r2,n,niici-(n.) 



Step 3. Estimate of 14. 
Indeed, ^4 satisfies the following equation: 

diY^ + Pld2Yi + iSldsY^ = -Y2d2B^ - Y^d^B^ - ys^a^o - ^353^0, 

y4(0,X2,X3) = 0. 
Then we have the following estimate: 

linii 

<Ci3<5f||yi||c2,.(f7.) + 11^2,^3,1411 



53S-^ 

(4.27) 



(4.28) 



(4.29) 



(4.30) 



Set Ci4 = maxjCg, Cii, C13}, then take e2 small enough such that Ci4e2 < 1. Now we 
set eo = min{ei, 62}. Then if < e < eo, A maps S to itself and is contraction in low order 
norm. Hence A has a unique fixed point, which will be the solution to ()4.8p . We have finished 
our proof. 

There are a few remarks in order. 

Remark 4.2. We remark that the vorticity of the back ground can be large. 
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Remark 4.3. The above iteration scheme can not be apphed to a 3-D nozzle with general 
sections. In general, the velocity field /32 and /Ss will lose one order derivative when integrating 
along the particle path. The delicate nonlinear coupling between the hyperbolic modes and 
elliptic modes {(^2,(^2) makes it extremely difficult to develop an effective iteration scheme 
in a general 3-D nozzle. How to effectively explore the good property of the quantity W = 
92/33 — d^h + kdi(32 — P2di(3s will be investigated in the forthcoming paper. 

5 A new formulation of the three dimensional incompressible 
Euler equations 



5.1 A new formulation for the incompressible Euler equations 

We apply the same idea to get a new formulation for the incompressible Euler equations. 
The original incompressible Euler equations take the following form: 



Define /3i 



Ui 



diUi + 92^2 + Ssus = 0, 

uidiui + U2d2Ui + usSsMi + dip = 0, 

UidiU2 + U2d2U2 + U'id^U2 + d2P = 0, 
UidiU2 + U2d2U2 + U3d3U2 + d2P = 0. 

1 



(5.1) 



— ,i = 2,3 and the Bernoulli's function B = —{ui + U2 + u^) + p- The inner 
ui 2 

product of u and the last three equations yields the Bernoulli's law: 



uidiB + U2d2B + usd^B = 0. 



(5.2) 



Multiplying the ffist equation in (15. ip by ^-^ , dividing the i-th equation in (15. ip by uf 

and adding them together for i = 2,3, 4, we obtain the following new system: 



dip2 + M3I32 - hd3P3 + —d2P = 0, 
diPs + 13282/33 - M2I32 + \d3P = 0, 



d2p2 + 93/33 



1 



(5.3) 



^dip = 0, 



diB + p2d2B + /3393S = 0. 
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Using the third equation, we rewrite the above system as fohows: 

1 



92/32 + 93/3; 



3 2 



dip = 0, 



5l/32 + /32 92/32 + /33 93/32 - ^dip + \d2P = 0, 



9i/33 + /32 92/33 + /33 93/33 - ^dip + \d:ip = 0, 



(5.4) 



9i5 + l32d2B + /3393S = 0. 

We also carry out a characteristic analysis to get an insight of the structure of incom- 
pressible Euler equations. 

We rewrite ()5.4p in a matrix form like the following 



Mi^iU + M292U + /W393U = 0. 



(5.5) 



Here U = (p, /32, ^3, Bf and 



-X 

ui 

'\l32 1 

■4/33 1 



Ml 



1 



,M2 



10 

\ (32 

/32 

/32 



.M3 



10 



/33 







/33 



/33 



Denote by A the root of det(AMi — ^2-^2 — ^3-^3)- A simple calculation shows that 



det(AMi - ^2^2 - 6/V^3 



1 



(A - /3 • ^f[X^ + (^2 + Ci)], then we find that it has one 



real root Aj. = /3 • ^ with multiplicity 2 and two conjugate complex roots A^ = ±iA/, where 



Xj = Y ^2 + ?3 • The corresponding left eigenvectors to Ar are the following 

1' = (0, 6 + W ■ 0, -6 - /32(/3 • 0, 0), 1^ = (0, 0, 0, 1). 
The corresponding left eigenvectors to A^ are 

1± = 1r ± ih = (-/3 • 6,6, 0) ± iiJ^Til 0, 0, 0). 



The differential operators corresponding to Ir and 1/ are {—{(3282 + (3383), 82,83,0) and 
(1,0,0,0). The action of these two differential operations will give us an elliptic system. 
We denote by = 82/32 + 83/33 and find that w and p satisfy the following elliptic system: 



-—081P + tZ7 = 0, 

ui 

8iw + 82{^82p) + 83{\d3p) - 82{^8ip) - 83{^8ip) 



(5.6) 



1 



+ /3393)(^9ip) + (92/32)' + {83/33)' + 282/3383/32 = 0. 
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Remark 5.1. The steady 3-D incompressible Euler system is always an hyperbolic- elliptic 
coupled system. Here our point of view is different from previous works Hy and jl29f . Indeed, 
in the works of f]^ and ]29^ . the authors regarded the pressure p as the Lagrange multiplier 
for the divergent free condition. 

5.2 The incompressible Euler equations with a constant Bernoulli's func- 
tion 



In the followmg, we assume that the Bernouhi's function is a constant and the speed |u| 

has a positive lower bound 6, then we have uf = As in the compressible case, 

the vorticity field will be parallel to the velocity field and W = 82/33 — 53/32 + fi-idij32 — /329i/33 
will take the place of 32^3 — d3U2. We define G = 1 + /3| + /3| and a new function K{p) 

of p satisfying K'(p) = — -. and then rewrite the equations satisfied by (32 and /33 as 

2{B -p) 

follows: 

G~\dil32 + P2d2/32 + - hdiK{p) + d2K{p) = 0, 

(5.7) 

G-\diP3 + p2d2h + M^h) - MiK{p) + d3K{p) = 0, 

Applying fi^di — 83 and — /32(?i + 82 to the above two equations and adding them together, 
then one can obtain the equation for W: 

BW + {82/32 + 83P3)W - G8iK{p)W - G-^WDG = 0. (5.8) 

Due to the first equation in (|5.3p . W satisfies the following equation: 

nW - G-^WDG = 0. (5.9) 

That is 

D — = 0. (5.10) 



G J 

It is easy to find G satisfies the following Riccati-type equation: 

DG - —^8ipG^ + —L^GBp = 0. (5.11) 

{B -p) [B- p) 

_i /3i 

Define three new variables Ki = G 2 ^Ki = — 7- = /3iKi,i = 2, 3. Here K2 and K3 behave as 

G2 

the sine function of "the flow angles". Define also two new functions I{p) = e s. 2t(B-t) 
rp I 

Q{p) = / ~77^ Tl{t)~'^dt, where p is a given positive reference pressure. Since the speed 

Jp 2(i3 -t) 

|u| has a lower bound 6, it is easy to verify that the above two integrals are convergent, so 

r(p) 1 

Up) and Qip) are well-defined. One should note that Up) satisfies — -— = — -. 

I VP) '^{B-p) 

It follows from (j5.1ip that Ki satisfies: 

T>Ki - — ^ rKiBp + — ^ rK-^8ip = 0. (5.12) 

2{B-p) ^ 2{B-p) 1 ^ ^ ^ 
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This equation, together with ()5.4I) yields: 
( 



1 



2{B-v) 



KiBp + diKi + d2K2 + 83X3 = 0, 



Di^2 
BK3 



2{B - 


P) 


1 




2[B - 


P) 


1 




2[B - 


P) 



K{Dp + 
K3'Dp + 



1 


2{B - 


P) 


1 




2[B - 


P) 


1 




2{B - 


P) 



K^^dip = 0. 
K^^d2P = 0, 
K^^d3P = 0. 



It follows from the definition of I{p) and Q{p) that 

1 



I{p) 
DiCi - 

DK2 - 
DK3- 



KiD/(p) + diKi + a2i^2 + = 0, 
-^i.,D/(,) + ^i.-a,/(,) = o. 

-^i.3D/(p) + ^K-a3/(p) = o. 



Setting Ai 



I{p) ' ' 



= 1,2,3, we obtain the following conservation laws: 

diAi + ^2^2 + ^3^3 = 0, 

di{Al) + 52(^1^2) + ^3(^1^3) + diQ{p) = 0, 
di{AiA2) + d2{Al) + 93(^2^3) + d2Q{p) = 0, 
di{AiA^3) + ^2(^2^3) + 53(^3) + d3Q{p) = 0. 
For the incompressible Euler equations, W satisfies the following equation: 



(Oi + P2d2 + Ms) 



w 



+ 



1 



G J ' 2{B-pf 
And the following conservation law also holds: 



[(l,/32,/33)-(Vpx Vi3)] = 0. 



di 



I{p,B) 



+ 82 



K2 



I{p,B) 



+ d3 



K3 



I{p,B] 



0. 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



Remark 5.2. As an application, one can apply the techniques developed in section^ to con- 
struct a smooth incompressible Euler flow in a rectangular cylinder, which satisfies the given 
incoming flow angles and the Bernoulli 's function at the inlet and the end pressure at the exit. 
Indeed, our background solution will be denoted by {pq,uq{x2,X3),Q,Q). The corresponding 
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Bernoulli's function will be Bo{x2,X3) = —Uq{x2,X3) +po- At the exit of the nozzle, we will 
prescribe the end pressure: 

p{l,X2,X3) = P0 + epeix2,X3). (5.18) 

Here it is required that pe{x2,X3) should satisfy the same compatibility conditions as (i4-4\ ) for 
Se{x2,X3) in the previous section. For the incoming flow angles and the Bernoulli's function 
and the slip boundary conditions, we should impose the same compatibility conditions as the 
compressible case. Then we can obtain the following existence and uniqueness results: 

Theorem 5.3. Given (/3^", /3|", G C^'^'iT^) satisfying the same compatibility condi- 

tions as (14. 2p and (14. 4p for the compressible case, there exists a positive small number eo, 
which depends on the background subsonic state {po,uo{x2,X3),0,0) and {/32"' , , B"^"^ , pe) , 
such that if < e < eo, then there exists a unique smooth subsonic flow {p,ui,U2,us)G 
C^'"(J]e) X {C'^'"{ne)f to dSH) satisfying boundary conditions (|ill . (f5lB and (jM]). More- 
over, the following estimate holds: 

||(ni,U2,n3) - (no(x2,2;3),0,0)||c2,a(n^) + |Ip - PolIc3."(Ce) ^ Ce. (5.19) 

Here C is a constant depending on {po,uo{x2,X3),0,0) and {f3^ , f3^ , B'^"' , pe) . 



6 Appendix 



In this appendix, we give the detailed calculations for (|3.4p - (|3.6p . 

Applying —f^^di + 53 and /32di — 82 to (|3.3p and adding them together, we obtain 

= {p^di - a3)(G-iD/32) - {Ml - 52)(G-iD/33) 

{Ml - damdi - d2)K{s) - {Mi - d2){Mi - d3)K{s) 

{Ml - 53)(G-^D/32) - {Ml - 92)(G-^D/33) - WdiK{s) 
J-WdiK{s). 



While 



J =G-'n{{Ml - 03)^2 - {Ml - 92)/33) 
3 

+ Y.^{Mi - d3)f3jdj(32 - {Ml - d2)Pjdj(3s 
i=i 

+ G-'Y.^A [(/^sSi - d3)dj - dj{Mi - d^m 

-[{Mi-d2)dj-dj{Mi-d2m 

+ {Ml - 53)G-^D/32 - {Ml - 52)G-iD/33 

1-1 



=G-'TiW + .h + J2 + h- 
Now we compute Jj, i = 1, 2, 3 respectively. 



(6.1) 



(6.2) 
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Ji = G-^ 



W{d2h + 53/33) - {Ml - 53)/3253/33 + {Ml - d^)M^h 



- {Ml - d2)M2Pz + {Ml - d2)M2h 

= G-\W {d2p2 + d^ip^) + Z]. 
Here Z = MiM^li2 - MiMsPs + MiM2li2 - MiM2k 

3 



J2 = ^/3, |[(/?35i - d^d^ - dj{Ml - 93)1/32 

j=l 

- [{Ml - d2)dj - dj{Mi - 92)1/33 



J3 = -G- 



G-^Z. 



{Ml - 93)GD/32 - (;329l - 92)GD^3 



3 r 



Pj{Mi - 93)/3,-D/32 - Pj{Mi - d2)PjBP3 



= -2G-^Y. 

i=i 

= -2G-2[t^(^2D^2 + ^3D^3) + J31]. 



Here 



J31 = MMi - d2)PzDh + Mhdi - d3)PzDh 

- MMl - 93)/32D^3 - MMl - 92)/32D^3 



= /32 



{Ml - 92)/33D/32 - (^29i - 92)/32D/33 



+ /33 



(^39l - 93);33D^2 - {Ml - d2)h'DPz 



h {Ml - d2)MiP2 - {Ml - d2)Mik 
+ k {Ml - d^)MiP2 - {Ml - di)Mik 



+ 



+ 



{Ml - d2)l5Md2h + Mzdzl32) - {Ml - 92)^2(^192^3 + Mzdzl5^) 



{Ml - d^)h{MM2 + /3i93/32) - {Ml - 93)^2 (/32/3392^3 + ^193/33) 



= GZ. 



(6.3) 



(6.4) 



(6.5) 



(6.6) 



Hence we have J3 = -2G'^\W{I32D(32 + (5^013^)] - 2G'^Z. Substitute these calculations 
into the above formula to get: 
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= J - WdiK{s) = Ji + J2 + J3 + G-^DW - WdiK{s) 
= G-^[W (82/32 + dsPs) + Z] + Q-^Z - 2G-^[W{^2D/32 + 

-2G-^Z + G-^DW -WdiK{s) 
= G-^DW + G-\d2P2 + dsPs - GdiK{s))W - G-'^WDG 

= G-^DW - G~^WDs - G-^WDG = pD{^). 

pG 

W 

This implies that {di + ^282 + MsK^) = 0. 

pG 
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